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I. INTRODUCTION 


In the past few years, improvements in the growth of epitaxial and 
bulk crystal forms of III-V compound semiconductor materials have led to 
many new devices. These devices include Gunn oscillators, light emitting 
diodes, infrared detectors, lasers, and special purpose transistors. The 
performance and lifetime of such devices depends sensitively on defects in the 
material. For example, the intensity of light emitting diodes is known to be 
strongly affected by dislocations that can be initiated by point defects. 

Most impurities in III-V compounds can be detected in very low con- 
centrations by chemical means. Vacancies, however, are difficult to detect. 

For this reason unexplained phenomena observed in these materials are often 
blamed on vacancies. It has been claimed that vacancies act as traps, recom- 
bination centers, and scattering centers. Since no method has existed to 
directly detect the vacancies, it has been difficult to confirm such assertions. 
Once low concentrations of vacancies can be detected, then their supposed 
effects can be confirmed; and, what is more important, crystal growth methods 
can be improved to eliminate them. 

The purpose of this work has been to develop experimental methods 
and related theory which will permit the measurement of low concentrations of 
vacancies and other defects in III-V compound semiconductors. Once the nature 
of these defects has been determined, this information can be incorporated into 
a transport theory for devices constructed from these materials , and experiments 
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conducted to test the theory. 

The vacancies and other defects in the Ill-V compounds are detected 

hy measurement of the nuclear magnetic resonance (NMR) line width. Most of 

thelll-V compoTUids have at least one isotope with a nuclear quadrupole moment. 

In a crystal with a cubic crystal field (characteristic of most III-V compounds) 

there is no quadrupole splitting of the Zeeman resonance line. However, a defect 

removes the cubic symmetry locally and causes splittings which result in a change 

of the NMR line width. This change can be used to detect the presence of vacancies 

l4 3 1 

in concentrations as low as 10 /cm . 

This report describes the progress to date on the project, as supported 
ly this research grant. Further reports will be made in the form of publications 
that are now being prepared for submission to scientific journals. Copies of 
these will be forwarded to NASA as they are completed. One of these reports has 
been submitted to the Review of Scientific Instruments for publication, and is 
included as Appendix A to this report. 



II. THEORETICAL ANALYSIS OF NUCLEAR RESONANCE LINE SHAPES 


The objective of this section is to compute the nuclear magnetic resonance 
(NMR) line broadening for the spins in a III-V compound semiconductor, in the 
large magnetic field limit. The spins interact normally via dipole-dipole 
and exchange interactions. In addition, charged point defects, such as vacan- 
cies or ionized donors and acceptors, break the cubic crystal symmetry and 
create a distribution of electric field gradients - and so q,uadrupole shifts - 
in their neighborhood. The electric field due to these charged defects shifts 
the ions around them from their eq.\jdlibri urn lattice sites into regions with 
finite electric field gradients. Such off-lattice site field gradients are 
large con^ared to those arising directly from the point charges . Since the 
electric field decreases only as the sq^uare of the distance r from each de- 
fect, each defect influences many ions. 

For low concentrations of defects where most ions experience small fields, 

—2 

the quadrupole shifts broaden the Zeeman lines. If the fields decrease as r 
the line shapes are Gaussian; if the fields arise from something like a dipolar 

_3 

charge distribution and fall off as fast as r the line shapes are Lorentzian. 
However, when the concentrations are higher, the shifts become large enough 
so that only the transitions between the levels identified by z-components 
of spin +1/2 and -1/2 remain in the original Zeeman line, and the intensity is 
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reduced. The vidth of this remaining line decreases because with these large 
quadrupole splittings some parts of the dipole-dipole and exchange interactions 
become non-secular. An experimental study of these effects will provide infor- 
mation about defect concentrations, and possibly the charge distributions about 
defects . 


Spin Hamiltonian 


The Hamiltonian for the system is : 


^2 " dJ " ‘ex 

The Zeeman Hamiltonian is : 


( 1 ) 


( 2 ) 


^ J. ■' -V 

\diere nucleus X has gyromagnetic ratio and spin . 

For like spins , and when Kq =0, the secTilar part of the dipole- 
dipole plus pseudo dipole-dipole Hamiltonian is: 




( 3 ) 


X/' 


Where the pseudo dipolar coupling constant arises from indirect interactions 

between the X and X nuclei coupled by p and d electrons . For unlike spins , a part 

-*■ — ► 

of contribute to the secular Hamiltonian. The changes when 

Ua ^ 0 will be discussed in detail later. 


The exchange Hamiltonian is : 

K. = tAA-Tr ’ 


2,£' 


ih) 


\Aiere is the exchange coupling constant that arises from interactions 

among the nuclei and s electrons.^ independent of the direction of the 

magnetic field, and thus produces isotropic changes in the line shape. 
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The secular part of the q^uadrupole Hamiltonian is: 

.2, 






i5) 


where is the quajdrupole moment of nucleus i , and the field gradient tensor 




is produced by the electric field 

.i 




at site 






k 

The indices i,J,k = x,y,z run over the principal cubic axes of the crystal. 

The piezoelectric strain terms and the direct field gradient terms have been 

J 

dropped from Q. becavise they are small compared to the field generated 

terms. The III-V compounds all have a T^ symmetry for which 



( 7 ) 

I I 

where R is the value of the only non-vanishing elements of the R tensor 

« S. Z Ji L I 

at site /,R=R =R =R .In the III-V compounds, R takes on 

3QTZ yzx zxy ’ 

only two values, one for Ill-atom sites, and the other for V-atom sites. 
Finally, from Eq. (6) we find: 

.i 


n.ec|'*'-n = i- 


( 8 ) 


Since most of the nuclei affected by an impurity lie outside the first 
few near neighbor shells, a reasonably accurate approximation for the electric 
field at displacement r from an in5)urity is: 

c - ix 
t - ^ 


( 9 ) 
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vhere e* is the effective charge of the impurity and e is the low freq^uency 
dielectric constant. Justification for such a simple form for the electric 
field in an insiolator is given hy Kubo.^ If, as we hope, experiments can be 
devised that provide detailed information about the charge distribution near 
an impurity, then a more accurate dielectric response formulation will be 
needed. However, useful insights can be gained from an analysis of this 
simple model. 

Since the. initial experiments have been conducted on GaAs, the remainder 
of this discussion will be specialized to the spin 3/2 case (the spin of 
Ga^^, Ga^^, and As'^^). The experimentally determined values^ of 

= 0.98 X 10^*^ cm"”^ and = 1.55 x 10^® cm**^ will be used. The natinral 
abundances, gyromagnetic ratios, and qmdrupole moments for the nuclei in GaAs 
are given in Table I. 


Nucleus 

Natural 
Abxindance 
(per cent) 

Gyromagnetic 

Ratio 

(kHz per gauss) 

Quadrupole 
Moment 
(cm ) 


Ga^9 

60.2 

1.0219 

- -2U 

.178 X 10 


Ga"^^ 

39.8 

I.298U 

-2k 

.112 X 10 


As'^5 

100 

0.T292 

-2U 

.3 X 10 




Table I. 



Calculation 

of the Number 

of Spins Influenced by 

a Defect 



Let us concentrate on one of the isotopes with Larmor frequency denoted 
hy V =Y^* magnitude of the shift of the Larmor frequency Av of the 
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m = I transitions from the states 3/2 to 1/2 and -3/2 to -1/2 of an ion 
z 

at displacement r relative to an imptirity is 


where 






eQe’*'R 

eh 




In order to compare our experiments with Rhoderick's^ experiments on the 
variation of KMR intensity with ionized impurity concentrations, or to compute 
the change of the line shape, the number of ions near an impurity whose quad- 
rupole frequency shift |Av[ is greater than some prescribed minimum 
must be found. To acconplish this we shall first calculate the volume near 

a given impurity where |Av|>Av^^. In order to take accomt of finite size 

-2 

effects and avoid infinities that arise from the fact that r falls off 

slowly, we shall actually compute the volvune of the sample for which 

M >Av , and eilso r<ro, where ro will be a characteristic dimension of the 

min 

sample. There will be two cases to treat. In the first instance Av is 
gTTfti 1 eno;agh so the volume is limited by ro> and in the second case Av is 
sufficiently lairge so the volume lies completely within the sample. 

This model for the finite size is a bit too simple, since it is correct 
only for an impurity located at the center of a spherical sample. A more 
accurate model would round off some of the abrupt changes this simple model 


introduces. We will find, however, that no observable physical phenomena 
are sensitive to this choice, so it is useless to deal with the complexities 
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of a more realistic model. The only purpose served "by the finite size is 

to provide a cutoff to eliminate an unphysical infinity. 

Another restriction must be placed on the lower limit of r, a<r, to 

account for the fact that the host atoms are separated from the inpurity by 

at least a distance a. For example, if the resonant nucleus is a type III 

nucleus, and the vacancy or inpurity is at a type V site, then a is the near 

neighbor distance a . If the defect sits on a type III site, then a is larger 

0 

(a = {8/3}^'**'^ ao). If the defect is in an interstitial site, then a will 
be less than ao. 

From Eq_. (lO), the condition that Av be constant is satisfied for 




where 6 is the angle between r and a. The broadening of the magnetic resonance 
lines and the decrease in intensity depends only on the magnitude of Av; 


positive and negative values contribute equally. Thus the shape of the volume 
with [ Av I >A v^^ ^ and a<r<ro is a figure of revolution about the direction a, 
of the shaded area shown in Figure 1. 
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For this case, the volume is given by the integral 


V 

= fff 


sin 0 ci0 


P'O 


4- il-IT 


- 


•CaS'KoJt'i/od 


Sin 0 d0 

0)S"’(ro*AV^') 


r^dr - 


(L 


r‘^4 


r^dr _ mr 


p^tt>S0/6V> 

Sin0d6 IrMr 

'^0 


(13) 




3 t) V- 6 oi. 

Next sii^pose that Av is large enough so that the volume lies within 

the sample. Then the volume over which the integral must be taken is a 

figure of revolution of the solid curve, for r>a, in Figure 1. Then the 

f%. r^Cose/A» 

integral is: ^ sinGdOj r^d r _ 

69 

If p is the density of nuclei of interest, e.g. Ga , then the number 

N (Av) with frequency shift greater than Av is 
s 

M^(av) = pXl(AV) 

0 


(15) 




where 


eQe*R 


and N is the number of relevant ions (Ga ) in the saii5>le. 


(16) 
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Calcination of the Line Shape Due to the Defect -Induced Quadrupole Interaction 

The distribution g(Av) of freq.uency shifts between Av and Av + d(Av) 
caused by a single impurity is related to N (Av) by the expression 


qCAV) -r ^ dNgCAv; 

^ dAV 


X - (a'^o ^ AV < A Vo 


( 17 ) 


= 1 ^ ^ (aV,/A.V)^^ - (aVoA'J^^i j AV„ ^ av < AV^ 


^ 0 , AV>AV^. 

Thus g(Av) is constant until Av reaches Avo, then decreases as (Av) 

2 2 2 

until Av reaches its maximum value (ro /a )Avo = (ct/a ), where g(Av) goes 
to zero. 

Next we must find the line shape and width to associate with a number 

of impurities distributed at random throughout the sample. For a broadening 

mechanism caused by dilute in5)urities, it is well known® that if the interaction 

between the impurity and the major constituent nuclei falls off faster than 
_2 

r then the line shape is Lorentzian. However, if the interaction falls 

_2 

off as r or slower, then the conditions for the "Central Limits Theorem" 

are satisfied and the line shape is Gaussian. Since the model we are investi- 

_2 

gating has the fields decreasing as r , we shall use the results of the 
Central Limits Theorem. 

2 

The second moment of the Gaussian Av"“ , according to the Central Limits 

2 

Theorem, eq.uals the second moment Av^ of the one impurity distribution times 
the number of impurities N^.. 


qChSiKAV) 

AV^ = ■ 


J g(AV)dA-^ 


~ K| (av [ i - ( AVq /aV>] 
^ 1 - (aV. /aV^T^ 




( 18 ) 
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where is the density of defects. As we stated earlier, Av does not 

depend on the size of the saii 5 >le, ro • 

Finally, the relation between the distribution of q.uadrupole freq.uency 

shifts G^(Av) (basically an Inhomogeneous broadening), and the NMR line shape 

G(v) must be established. The quadrupole frequency shift distribution is a 

2 

Ga\;issian function with variance Av : 




(19) 


For the spin 3/2 case, the magnetic resonance line shape is found through 


the expression 

J>0 


(9i(v) = 3 ^^ C>?-AV AV)4- ’*‘.3^3^ AV^ ^ 

^ ( 20 ) 


where G . ... (v -Av, Av) is the line shape function of the 3/2-»-l/2 transl- 

tion caused by other broadening mechanisms (as modified by the quadrupole 
shifts). The first argment in the functional dependence of G 3 ^^ is the 
center frequency, and the second is intended to represent the effect the shift 
may cause on the line shape. The shape due to other interactions can be modi- 
fied by the shifts, e.g, terms in the dipolar Hamiltonian that are secular in 
the abscence of shifts will become non-secular when the shifts are large 
enough, thus modifying the width of the dipolar broadened line. G_i, i, and 
G_^^^_,^^are the corresponding shape functions for those transitions. The 
factors 0,3, O.U, 0,3 multiplying the shape fimctions are weighting factors 
which arise from the transition matrix elements for the three transitions. 

We shall have to treat the more complex case later, but for now examine 
the simple situation corresponding to low defect concentrations. 
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2 I )Z 

I ~ (^1 I 

‘i)2 ' "i,2 


C ,(v + iV/V)=:4 , (V^o) 

t r 2- 


J.) 2 


I =Q. . • 

2-j 2 1)2 J ) Z 

Then Eq. (20), in view of the approximation in Eq. (2l), "becomes 

OO 

C:r(y) ^ 0,(^c/^)) G>^()^-^6 '\>) {a^} . 

•^OO 

If G (v) is a Gaussian with variance Oj then 
0 

, 0 .^ -P.<T,^ 

~ r:r::-f^i. ~2 Vi C f- C 




^zif % 


( 21 ) 


(22) 


(23) 


where V is centered at the Larmor frequency. Thus the variance of the 
composite line is 



(o^^+ 0.4, AvO 



(2U) 


The contribution from the quadrupolo interaction can also be found 

from a standard "Van "Vleck calculation that yields an answer involving a lattice 

sum.^ If one uses the electric field-induced gradients, as given in Eq. (l6), 

2 

the angular dependence of CJ^ is reproduced exactly. If, in addition, the radial 
lattice sum is replaced with an integral, the same numerical result, Eq. (24), is 


obtained. 



III. EXPERIMENTAL EQUIPMENT 


1. Pulsed NMR 

The descriptions of equipment in this section may be clarified by 
reference to Fig\ire 2, a block diagram of the pulsed NMR system. Much 
of this system was purchased with funds provided by a National Science 
Foundation Departmental Development Grant. However, graduate student support 
and faculty summer salaries provided by this grant were indispensable to the 
construction of the system. 

a. Pulse programmer 

A programmable piilse generator was designed and built in this 
laboratory to control the sequences of pulses of radio frequency (rf) magnetic 
field that are applied to the sample. A number of pulses and delays can be 
called for serially in time. Each pulse can have one of four different rf phases 
in the gating and phasing circuits. The rf can be supplied either as normal 
high intensity pulses, or as low and variable intensity pulses for "rotating frame" 
experiments. Both the pvilse d\iration and the delay time can be controlled either 
digitally or by analog circuits. 

Although the pulse programmer is limited to l6 separate instructions, 
two capabilities extend its effective memory. Many pulse sequences are repeti- 
tions of a basic short sequence. The programmer has two separate counters that 
permit "looping" of a particular sequence a number of times. In addition, the 
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Figure 2. ELock Diagram of Pulsed NMR Apparatus 
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pulse programmer memory can be loaded by the minicomputer that controls the signal 
averaging process. Thus the programer memory is effectively extended to the 
memory available in the minicomputer, if there is enough time (approximately 
50 microseconds) in a given sequence to reload betveen pulses. 

A more complete description of the programmer is being prepared for 
publication (see Appendix B). 

b. RF Phasing and Gating 

This portion of the pulsed NMR system provides piilses of rf in 
the sequence specified by the pulse programmer. The rf source is a General Radio 
frequency synthesizer with excellent frequency, amplitude, and phase stability. 

The pulses may have any one of four phases with respect to the master oscillator; 

0°, 90°, 180°, 270°. 

Actual pulse gating and amplification are done at a fixed frequency, 
and the result is mixed with a variable frequency to provide the NMR stimulating 
pulse. This procedure eliminates much of the tuning usually required to change 
NMR frequencies. 

c* Power Amplifier and Sample Probe 

A combination of commercial power amplifiers and circuits made in 
this laboratory is used to apply the rf pulses to the sample* The main power 
amplifier is a NMR Specialties P-103R> capable of applying up to five kilowatt 
rf pulses. 

ft 

The "single coil" probes that are used combine suggestions by Waugh^ 

9 

and Lowe . The rf magnetic fields produced are about 75 gauss, which makes the time 
for a "7ry2 pulse" for about 3-3 microseconds. Recovery time of the system 

after a pulse (until the NMR signal may be seen) is on the order of 3 to 5 micro- 
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seconds at 20 megahertz. For Ga in GaAs, the NKR signal decay time (T^) is 
about 150 microseconds. Thus our pulse amplitude and recovery time have been 
more than adeg^uate for these experiments . 

d. NMR Signal Detection and Averaging 

The transient response of the nuclear magnetization after the 

stimulating rf pulse is amplified and detected by a combination of commercial 

8 9 

apparatus, and circuits like those described by Waugh and Love. The signal 
may then be viewed on an oscilloscope screen or undergo analog-to -digital 
conversion. It is this latter process that provides the large signal -to-noise 
ratio required for our experiment . 

We have constructed a digital processor and computer interface to 
provide the following signal averaging process: (l) analog-to-digital conversion, 

(2) temporary storage of data in a high speed buffer memory, (3) transfer of 
data to the memoiy of a small computer, (U) averaging and simple manipulation 
of data by the computer, and (5) transfer of data to IBM-compatible magnetic 
tape for use on a large scale computer. 

The analog-to-digital converter samples the NMR signal at a maximum 
rate of once every 200 nanoseconds, and converts the amplitude to an eight-bit 
binary result. This result must be temporarily stored, since the data rate is 
too high for present-day con^uters. Timing and data transfer are accomplished 
with digital circuitry designed and constructed entirely in this laboratory. 

This circuitry will be described in detail in a paper now in preparation (see 
Appendix B) . 
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The small computer is a Digital Equipment Corporation PDP-8/E. It 
has internal Interfacing to allow display of memory contents on an oscilloscope, 
and to record data on a cartridge -style magnetic tape recorder. Simple addition 
of incoming data to the previous data stored in the computer improves the signal- 
to-noise ratio, since the signal amplitude adds up more rapidly than the noise. 

Memory locations currently available in our computer allow us to sample 
up to 102^ points per transient signal. This provides a very detailed picture 
of the NMR signal. This signal can be Fourier transformed to show the usual 
frequency domain line shape, or manipiUated directly to find certain average 
values, such as the second moment. 

Another feature available in our interface allows analog-to-digital 
conversion to be done at selected times during the decay. This aids in analyzing 
the complicated structiire of the signal resulting from a pulse sequence. 

Finally, we have constructed a simple interface that permits us to 
use an IBM-compatible magnetic tape deck for recording the data from our computer. 
We are then able to transfer the data to the William and Mary IBM 36 O /Model 50 
conputer, and to take advantage of the library facilities and peripheral 
equipment of a full scale computer center. This allows, for example, display 
of the data by a CalComp plotter-see Figures 3 €Uid 4, A more complete description 
of the interface has been given in a paper submitted to the Review of Scientific 
Instriiments (Appendix A). 

2. CW MMR 

We have available an KMR spectrometer built by D. R. Torgeson,^*^ and a 
Varian V-it230B crossed-coil rf probe . This combination has an excellent signal -to- 
noise ratio, and samples and frequencies are easily changed. We also have a 
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Northern Scientific model NS -561 Signal Averager. This unit has also been 
interfaced to the small computer, and thus may have its data transferred to 
the IBM system, 

3. Masnet 

Fields of up to i 9 .lt kilogauss are provided by a Varian model V-3900 
electromagnet and power supply. The field is stabilized with a Varian Mark I 
Fieldial Hall effect device. This system was chosen for its combination of 
strength, homogeneity, and stability of field, all of which are essential to 
the experiments proposed here. 



IV. ANALYSIS OF DATA 
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In this section, the results of Ga NMR studies on two different 
samples of GaAs will be discussed. This represents the data that had been 
analyzed by the completion date of this grant j however, far more data have 
been taken and analyzed during this past summer. Those data will be briefly 
described at the conclusion of this section. 

Figures 3 and h are computer plots of typical transient NMR signals 
from the two samples. This signal, usually called the free induction decay 
(fid), is actually composed of 102k separate data points that have been 
accum\ilated by the signal averaging system discussed in Section III. The 
plotter has interpolated between points , making the signal appear continuous . 

69 

A pure, undamaged sample of GaAs gives the Ga FID shown in Figure 3- 
The signal averaging process is actually initiated before the rf pulse, giving 
rise to the horizontal baseline to the left of t=0. At t=0, the rf pulse is 
applied, and the FID begins after the pulse. There is a short dead time (about 
10 psec . ) , and some noise generated by the piezoelectric response of the crystal 
at the very beginning of the FID. The remainder of the signal is nearly Gaussian in 
shape, but with a pronounced oscillation - seen most clearly in the negative-going 
curve around t=U80 psec. 

Figure U is the FID of Ga^^ in a GaAs sang>le that was heated to 800°C 
to create vacancies, then quenched rapidly to room temperature to trap them. 
According to Pearson, Potts, and Macres,^^ the number of vacancies present should 
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FID AMPLITUDE 



Figure 4* FID, Ga ^As in Damaged Sample 
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be about 5 x 10"^° This cvirve appears more exponential in character than 

Figure 3» and the oscillation is smaller. 

Similar curves were obtained for these samples as a function of the 
orientation of the sample in the magnetic field. The data were transferred to 
magnetic tape, then processed in the William and Mary IBM 360 computer. 

lieast— scjuares best fits were obtained for these data, assuming a 
modified form of the function given in Eq. (23). The form used was 

Gft) = c(,6 +.‘i ■ (25) 


The essential difference is that the experimental results are expressed in the 
time domain, while Eq. (23) is in the frequency domain. The curve fitting is 
quite successful, as seen in Figures 5 and 6. 

In order to compare the theory to these results, Eqs. (l8) and (2U) 
are used. The second moments of the experimental results are found from Wlf.O 

Figures T and 8 display the experimental second moments of the two 
crystals as a function of crystal orientation in the magnetic field. Also shown 
are the Van Vleck calculations of the magnetic dipole-dipole second moment, 
and a best fit curve of the functional form 


= A + B fe) = A + 6 & “s'® - 1 • 


(26) 


This function is similar to Eq. (l8), and should have the correct angular 
variation for all of the NMR line broadening mechanisms present in this material. 
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Line broadening mechanisms 

2 

As discussed by Sundfors, there are three line broadening mechanisms 
important to the data taken in these samples. They are: magnetic dipole-diple , 

electric quadrupole with impurity -caused field gradients, and nuclear pseudo- 
exchange interactions. In the undamaged crystal, our calculation of dipole- 
dipole and exchange effects are in agreement with Sundfors. Further, the 
exchange effects are not different in the damaged crystal. The remaining questions 
are thus: how does the dip ole -dipole contribution vary with defect concentration; 

and more importantly, how sensitive is the quadrupolar broadening to the defect 
concentration? 

The dipole-dipole interaction varies with defect concentration due to 
the first-order quadrupole frequency shift of the NMR 3/2-^l/2, -l/2“> -3/2 
transitions. This shift makes part of the dipole-dipole hamiltonian non-secular, 
reducing the line broadening and causing an apparent narrowing of the resonance 
line. This effect is a small change - 10% or 20% - and we have discussed it in 
an earlier report on these experiments.^ 

Most of the change in second moments between Figures 7 and 8 is due to 
quadrnipole mechanism proposed in Section II. The enormous change due to charged 
defects verifies our hypothesis that NMR is a sensitive tool for detection of 
these defects. 

To be quantitative in this analysis, one must change the direction 

cosinos expressed in Eq. (l8) to the angle of orientation of a sample in the 

2 Z "2. 

magnetic field. First, using the general expression — JL 

J 0 

one can show 
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These samples were rotated about the ^10^ axis (i.e. the axis 
is held perpendicular to the constant magnetic field.)- Defining the angle 
between a^|To)axis and the magnetic field as 0 , Eq. (l8) becomes 



(18') 


where f(6) is the function defined in Eq. (26). The quadrupole contributions 
to the second moment have the same functional form as the ang\ilar dependent part 
of the dipole-dipole contributions. For example, the Van Vleck calciilation of 
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dipole-dipole broadening for Ga As gives 


<AUi > = .29<? -^ / 53 / 'S') (k H/) . 


In order to find the quadrupole contribution to the composite 


(2T) 


second moment <X! 


, one must subtract all other contributions from the composite 
moment, as shown by Eq. (2U). Adding Sundfors ’ value of the isotropic exchange 
coupling, ^ “ , 2 ,( 0 ^ ^ dipole term given by Eq. (27), 

the sum of terms other than quadrupolar becomes 

=.5C‘f4-U3-f(e). 

The composite second moments for the -undamaged crystal, 
the damaged crystal, given by the best fit values for the 

parameters A and B in Eq. (26). The cxirves shown in Figures 7 and 8 are 


and 


4- iZZ f(e) (undamaged) , 


(29) 


and 


( ff""! = Z.48 4 4.30 f(e) (damaged) . 

^ fri 


} 
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Using Eqs. (24), (28), and (29), one finds 

^^2 ^ ,Q 

and (30) 

“ 3,zo + 4^,^5f(e ) . 

If no other hroadening mechanisms are present, the coefficient of f(G) 
in Eq. (30) arises from quadrupolar broadening. Two aspects of Eq. (30) 
make this attribution suspicious - the existence of an isotropic contribution. 


and the negative coefficient of f(©) for the undamaged crystal. The latter 
implies that the dipolar and isotropic exchange interactions give contributions 
larger than the experimentally measured moments. . More recent data have shown that 
a pseudo— dipolar exchange interaction is present. This will be reported at a 
later date. For now, for order-of-magnitude purposes, we will assume the 
coefficient 4.95 in the damaged sample is due entirely to quadrupolar broadening. 

The angular dependent part of gives us in Eq. (l8*)- 

Using the values (X ~ (assuming an As defect), =• IZ.S ^ 

^ \0 A — 

R= lO Cwi"'* , and Q from Table I, we calculate lO . 

This value is an order of magnitude smaller than that predicted by 

, 11 / 3 N 

extrapolation from the data of Pearson, e^.^. I = 5 tO )• 

Considering the crude analysis here, this discrepancy is not at all disturbing. 

The important point is that we have easily detected a defect concentration that 
is lower by a factor of at least 100 from previous measurements. 

More Recent Data and Analysis 

During this past summer, considerable additional data were taken. 

Ga^^As, Ga"^^As, and GaAs"^^ resonances were observed as a function of crystal 
orientation; for an undamaged, pvire crystal, a crystal doped with Si (carrier 
concentration 4x 10^^ cm~^), and pure crystals heated to 500°C, 600°C, and 
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700°C and then q.uenched at room temperature. Signals were observed carefully 

immediately after quenching in order to observe any annealing effects. 

Using 1 three resonances, one can calculate the contributions of 

all broadening mechanisms . The quadrupole broadening thus calculated gives 

impurity concentrations that generally agrees in magnitude with the data of 
11 

Pearson et al> , but at even lower concentrations than cited here* 

As the analysis of these data is completed, and the results prepared 
for publication, N.A^S.A* will be sent copies as addenda to this report. 
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Abstract 


A sin5>le and inexpensive computer interface for vrite-only parallel 
data transfers from a Digital Equipment PDP-8/E miniconputer equipped with a 
KA8-E Positive I/O Bus to some peripheral device is described. The data trans- 
fer is under program control and is thus versatile enough to transfer data to 
a nvunber of devices. The application discussed here is a data transfer to 
magnetic tape xising a Cipher Model 100 incremental tape drive, for later 
processing by a large conputer. 



It is often advantageous to collect data using a laboratory minicon 5 >uter 
and then to transfer the data to a larger con^juter for analysis. In this laboratory, 
signal averaged pxilsed NMR data are recorded on IBM con 5 )atihle 9 track magnetic 
tape by our Digital Equipment Corp- (D.E.C.) PDP-8/E, using the interface described 
here. The interface allows parallel data transfer under program control and is 
able to coviple the PDP-8/E through the Positive I/O Bus KA8-E to many devices 
merely by changing the output cable. 

We record data teirq)orarily at the end of each run using a Tennecon^ TP-1375 
cartridge recorder, and transfer it as a group to IBM compatible tape at a later 
time. Thus the tape drive is only occasionally needed and can be shared with other 
research groiJ^s, In addition, punched paper tape that is used by many laboratory 
instruments can be read into the PDP-8/E through the associated teletype, converted 
to the appropriate format, and transferred to the magnetic tape. 

The circuit design, shown in Figure 1, requires eleven integrated circuits. 
Three MC3001 are used to buffer the data, preventing possible loading and cable noise 
problems. As designed, data are always present on the peripheral inputs, but if 
necessary the twelve AND gates coizld be enabled by a buffered input/output pulse 

i 

(BIOP) from the KA8-E. 

An eight input HAND gate (MC3015) detects the correct device code (as 
designed, 2Ua) specified by the buffered memory bits 3 throxi^ 8 (BMB 03-08). The 
NAND gate output enables the instruction decoder, a one-of-ei^t decimal decoder 
(MCUoUS). This design uses only the instructions binary coded 1 throu^ 4 by 
BMB 09-11. Coincidence of an instruction decoder output and a BIOP pulse on the 
input of one of four AND gates (MC3001) executes the instruction. 


1 



Two instructions may assert the KA8-E SKIP L signal (in turn causing the 
PHP-8/E to skip the next program step). The peripheral BUSY signal (af^er inversion) 
inhibits BIOP 1 from asserting SKIP L if the program interrogates the BUSY status 
(instruction l) while the peripheral is busy. In addition, if the peripheral 
detects an error, a positive level on one of four inputs is inverted and passes 
through the MD gate (MC3011) . This allows BIOP 2 to assert SKIP L when the 
program interrogates the error condition (instruction 2). Four lamps and drivers 
(MC 858 } are used to identify the error. Open collector NMD gates (N 888 I) assert 
the SKIP L line . All inverters are in MC3008 hex inverter packages . 

Instruction 3 catises a STEPA^RITE action by the tape drive. Finally, when 
the data transfer is complete, instruction k causes the tape drive to write an inter- 
record gap ( IRG) . 

The External Bus cabling plugs into a D.E.C, H803 connector block. This 
interface \ises half of a D.E.C. W979 module. The output cabling consists of twisted 
pairs anchored to another module that plugs in an adjacent slot in the block. 

The recording program depends on the type of perijhered.. Most parallel 
data transfers can be done with a combination of the fovir instructions given here. 

For example, this same interface is being used to transfer data to a Digi-Data Model 
1700 synchronous tape drive and formatter, and directly to an IBM -360 equipped with 
a special IBM 2972 Scientific Interface Control unit. 

Copies of the software, both the PDP-8/E program and the Fortran program for 
reading the magnetic tape, are available from the authors. 



KA8-E CIPHER 

EXTERNAL MODEL 

BUS 100 



Fig. 1. Interface circuit diagram and interconnections. Cipher tape 
drive error signals are (l) broken tape, (2) tape end, (3) file protect, and 
(I*) parity. Not shovm are several inputs to the tape drive that must be 
gro\mded for write-only operation. 
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ABSTRACT 

Piilsed NMR was used to detect charged defect 

l4 -3 

concentrations as small as 2x10 cm . An undoped type n 

single crystal of GaAs was cut into several samples, and 

NMR second moments were obtained as a function of magnetic 

field orientation for Ga"^^, Ga^^, and As”^^. The orientationally 
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dependent part of the Ga second moment is one half that 
expected from the dipolar interaction. An explanation is given 
based on interference between the dipolar and the negative 
pseudodipolar interactions. There is a contribution to the 
second moment due to the electric quadrupolar interaction 
between nuclei and the electric field gradient (efg) associated 


with point charge like defects, and it is proportional to 
the defect concentration. Because there are three isotopes 
we can separately identify the quadrupolar and pseudodipolar 
second moments. In order to introduce crystalline defects, 
samples were held at a constant temperature (500°C, 550°C, 

600°C, 700°C) in evacuated tubes and quenched to room 
temperature. The increase in second moment is due to the 
increased quadrupolar contribution, and it establishes the 
defect density for each damaged .sample. The pseudodipolar 
interaction is observed to be independent of damage. Cto 
data indicate that the defects are located on the As nuclear 
sites. Ko annealing was observed. 

*Supported in part by NASA grants NGR-U7-OO6-O5O and 
NGL-47-006-055 and NSF grant GH-4l082. 

"^Based on a dissertation submitted by one of the authors (R.K.H. 
the Department of Physics, College of William and Mary, 
in partial fulfillment of the requirements for the degree 


of Doctor of Philosophy- 



I. INTRODUCTION 

In III-V compoimds, the broadening of the nuclear magnetic resonance 

(NMR) line is due to four rigid lattice interactions. In particular, the 

NMR second moments of the three isotopic- species in GaAs (Ga , Ga , and As ) 

contain contrihutions from all four of these mechanisms. 

The dipolar interaction is well understood, and its contribution to 
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the UMR second moment may be calculated theoretically - In 1957 Shulman, et al. ’ 
reported linewidths in un doped powder samples of GaAs larger than those expected 
from the dipolar interaction alone. The additional broadening was attributed 
to the indirect exchange interaction between unlike spin systems. This 

nuclear exchange interaction involves an indirect nuclear spin coupling via 
the hyperfine interaction between the electronic and nuclear spins. Its magnitude 
depends mainly upon the s character of the electronic wave function at the 
nuclear positions. More recently, Sundfors”^ has observed linewidths in several 
III-V compo\mds using both NMR and nuclear acoustic resonance (NAR) techniques. 

He identified the exchange contribution to the second moments and, using the 
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theory of Anderson , the exchange coupling constants . 

Another electron coupled interaction is the pseudodipolar interaction . 
The role of this mechanism in the broadening of NMR lines in III-V compounds is 
not well established, but a large negative pseudodipolar effect (canceling some 



dipolar broadening) has .been observed in the resonance in InP . The 

pseudodipolar interaction will be required to explain the data reported here. 

In undamaged samples both the Ga and Ga NMR second moments are smaller than 
those predicted by the dipolar interaction alone. In this case, the sign of 
the pseudodipolar coupling constant may he unambiguously determined, 

GaAs has a zinc blende structure. Thus, the efg at nuclear sites is zero 

unless there are lattice defects. Randomly located lattice defects cause a random 
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efg at nuclear sites and broaden the NMR line. E. H, Rhoderick observed the 
effects of suhstitutional impurities upon NMR line intensities. The intensity 
loss with increasing impurity concentration is due to quadrupolar broadening 
caused by the efg associated with the impurities. It was found that the efg 
associated with ionized donor or acceptor impurities is larger than that 
eissociated with strains produced by neutral impurities. The explanation of the 
sensitivity of the NMR line to the ionized impurities^^ required anomalously 
large antishielding factors on the order of 1000. 

Gill and Bloembergen^^ reported a direct measurement of the coupling 
constant relating the components of an efg tensor to an applied homogeneous 
electric field. The efg induced at nuclear sites is caused partly by a distortion 
of the valence orbitals and partly by a relative displacement of the Ga and As 
sublattices. Taking both effects into account, the coupling constant was 
explained using antishielding factors of 2h and 30 for Ga and As, respectively. 
Since the induced efg is proportional to the electric field, its magnitude 
decreases as r in the field of a point charge. In contrast, the direct efg 
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caused by a point charge decreases as r • Thus the quadrupolar 
interaction is dominated by the induced efg. Using this induced quadrupolar 
interaction it is not necessary to introduce large anti shielding factors to 
explain Rhoderick’s data* 

Sundfors^ identified the electric q^uadrupolar contribution to the NAR 
linewidths in samples of known carrier concentration* Comparing his measured 
linewidths to those predicted by a simplified theory of the induced 

q.uadrupole interaction linevidth, in which the number of charge centers is set equal to 
the carrier concentration, he obtained order of magnitude agreement with the 
measurements of Gill and Bloembergen* 

In this research we demonstrate that WMR linewidths can be used to 
detect low concentrations of charged defects in single crystals of GaAs* We 
have done this with careful measurements of NMR second moments as a function of 
magnetic field orientation for all three isotopic species* We identify individual 
contributions from the four broadening interactions; and, using an improved theory, 
we show that the first order quadrupolar contribution is proportional to the 
charged defect concentration* Thus the identification of the quadrupolar 
contribution to the second moment is a measurement of this concentration* 

Defect concentrations in several samples were altered by subjecting the 
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samples to thermal damage* The process was similar to that of Potts and Pearson , 
but our damage temperatures were far lover, producing fewer defects (on the order 
of 5x10^^ cm"^, compared to their 10^® cm"^.) Using vapor pressure data we estimate 
the As vacancy concentration caused by As sublimation at elevated temperatures* 

In this estimate we assume that As molecules in the vapor phase are frozen 



out of the lattice by the thermal quench. The vacancy concentration deduced 
from the vapor pressure data is more than an order of magnitude larger than the 
concentrations projected from Potts and Pearson^ s experiments. Our NMP measure- 
ments of the defect concentrations lie about mid-way between the two projections . 



II . THEORY 


Consider "two different nuclear— spin systems with nuclei having spin 
angular momentum *KX and*lfvS , and non-overlapping resonances. The resonance 
experiment is performed on the spins I. We assume that the resonance is 
broadened by the following four interactions: dipolar, pseudodipolar , exchange, 

and quadrupolar. The second moment of the resonance line is calculated from the 
expression given by van Vleck^ 

^ TrLW,iI/Tr[t*], (1) 

where is the system hamil’bonian , and I is the total spin operator proportional 
to the observed magnetization. 

Following the notation of Sundfors , we write the system haimltonian as 






( 2 ) 


where and *5^8 are the Zeeman energies for the I and S spins, respectively. We 
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may combine the truncated dipolar and pseudodipolar interaction terms as 


X:* X'- Pxe.,) ^ , 


(3a) 


Loi 0 5. spPX»>i) nr , 


(3b) 



(3c) 


/ 

where B. . is the pseudodipolar coupling constant for like spins I, B is 
pseudodipolar coupling constant for like spins S, is the pseudodipolar 

coupling constant for unlike spins I and S, and r.^ is the displacement vector 
from the i^^ nuclear site. The second Legendre polynomial is written as 



where Q is the angle between r. . and the Zeeman field. The exchange terms may 
vij = 

be written 



(3d) 


(3e) 


(3f) 


yliere A., is the exchange coupling constant for like spins I, 

exchange coupling constant for like spins S, and is the exchange coupling 

constant for unlike spins I and S, The quadrupolar term can be written 


vtiere V . is the efg at nuclear site J, and A=eQ [^l(2I-l)] , where is the 

zzj 

electric quadrupole monient of the spins !• 

By combining Eq^s* (1)» ond (3) one obtains the expression for the 

NMR second moment o ^ 2 ^ 

<4H‘> = 31(1+1)^*^^ +iS(5t02l^Aij 



The lattice sums in the dipolar terms of Eq. (4) may "b^ separated 
into tvo parts, one of which contains no pseudodipolar coupling constants. 
Consider just the dipolar sum for like spin broadening. This may be written 


The first siim on the right side of Eq. (5) contains the pseudodipolar contribution 
to the second moment and the interference term between the dipolar and pseudo- 
dipolar interactions. The second sum is just the dipolar contribution in the 
presence of zero pseudodipolar interactions. This sum may be performed, and for 
the zinc blende lattice it has the form^ 


n ; C SIX') = 5 + b(v«3+v^?+ m;)] ^ 


( 6 ) 
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where n . n , anti n are the direction cosines of the Zeeman field in the 
x’ y z 

crystalline coordinate system. When the sum is performed over nuclear sites on 

the fee sublattice containing the spins I, Sundfors computes a=256 and 

a where a is the lattice constant. For the non-resonant sublattice 
o * o 
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he computes a *=2491 Q-nd b*=-2390 9-^ 

If the magnetic field is rotated in the (llO) plane, and 0 is defined 
as the angle between the [110] direction and the Zeeman field, then 


mJ + riy + T^i = 1“ *^(0)^ 

where f(6) = 2 cos^ - j cos^Q . Combining Eqs . (6) and (t) we may write the 
dipolar sum in terms of f (9) , 


( 7 ) 




( 8 .) 



where X = (a+b)/8, and Y - -*b/8. The dipolar contrihution to the NMR second 
moment thus contains an isotropic part, X, and a part dependent upon the magnetic 
field direction, Y. 

We consider now the first sum on the right hand side of Eq* (5)* 
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electron coupled pseudodipolar interaction is extremely short range , and we 
assume that only nearest neighbors interact* Thus the first sum in Eq* (5) 
vanishes because it does not contain the nearest neighbor shell. However, the 
corresponding tenn in the expression for 'unlike dipolar broadening does contain 
the nearest neighbor shell. The coefficients B. . are independent of 0 and may 
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be factored out of the sum. The pseudodipolar terms become 


f(e)^ 


(9) 


where is the nearest neighbor pseudodipolar coupling constant. With this 

nearest neighbor approximation, the pseudodipolar terms have the same orientation 

dependence as the angular part of the dipolar contribution to the second moment, 

but there is no isotropic pseudodipolar contribution. 

In order to evaluate the quadrupolar contribution to the second moment 

we require a specific model for the efg tensor. The efg associated with crystalline 
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strain is much smaller than that caused by charged impurities ’ and may be 
neglected when both are present. Thus we assume that the efg is due entirely to 
the electric field of charged crystalline defects. The general tensor relationship 
between the components V. of the efg tensor and the electric field con®)onents E 

X J 


is given by' 


.11, U 



( 10 ) 


^i.1k 

system Eq. 


is a third rank, tensor. In the zinc blende 
(lO) may he written^^ 


crystalline 


coordinate 




/o £, £^\ 

Eli O Ex 
\ Ey £Tx O j ^ 


(11) 


where is the magnitude of the non-vanishing tensor components. 

The component of V in the secular part of the quadrupole hamiltonian 
is n^V'Hj where n is the unit vector in the direction of the Zeeman field. We find 

n- — 2Ri‘i(T\Y'^E'k + 

Since most of the nuclei affected by a given impurity lie outside the first few 
near neighbor shells, we use the static macroscopic dielectric constant €. to 
approximate the electronic screening. Thus the electric field at a displacement 
r from a single impiLrity of effective charge e* is 

E = /er^ (13) 

Combining Eqs. (!•), (12), and (13), we obtain an expression for the quadrupolar 



contribution to the second moment for a single impurity 




Due to the lattice symmetry of GaAs, summations over terms 
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vanish, and summations over terms x., y., and z. are equal. Using these and the 

*3 t3 J 

properties of direction cosines ve may rewrite Eq. {ih) 


<aw% = I (-f^T . 


( 15 ) 


We assume that the charged defects occur randomly on one sublattice. 

The efg at a nuclear site is the sum of contributions from each of the charged 
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defects. From the central limits theorem, the broadening due to impurities 
is Nj times that due to a single impurity, and the quadrupolar second moment 


expression becomes 


<4 f [SKl+0-3] (t^) ' , 


(l6) 


where the ratio of defect density to nuclear density p^/p has replaced the ratio 
Nj/K. This contribution has the same angular dependence as the pseudodipolar 
second moment. Both vanish when 0 =1T/2, where the Zeeman field is aligned with 
the [OOl] direction. 

We have calculated the lattice sum over r.^ in Eq. (l6). The discrete 

SUM was p6rf‘onned over a cube with an edge size of l6a^* An integraJ. was perfomed 

XT • ♦ * — U 

over the remainder of the sample . The result of this calculation is 33*75 

for impurities located on the same suhlattice as the nuclear spins I, and 63-31 

a for the impurities located on the other suhlattice* 



Ill . EXPERIMENT 


The experiment was performed on several type n single crystals of 
GaAs, sliced from a single imdoped, hoat grown ingot, purchased from the Monsanto 
company. The properties of the ingot were supplied hy Mansanto as follows, resistivity, 
1. 3^x10^ ohm— cm; carrier concentration, 2.4x10^^ cm electron mobility, 2.7x10 
cm^volt"^sec“^; and etch pit density, 4.2 x 10^ cm”^. A 50 micron diamond 
saw was used to cut the ingot into several samples of nearly equal size, about 
0.5cm by 0.5cm by 1.0 cm. A type n silicon doped sample purchased from 
Electronic Materials Corporation was also studied. Its properties as 
measured by the manufacturer, were: resistivity, .03 ohm-cm; carrier 
concentration, 4xl0^^cm electron mobility, 5x10^ cm volt sec ; and etch pit 
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density, 8x10 cm . ' 

The free induction decay (FID) following a single 90° pulse was observed 
for the three isotopes in each of the samples as a function of magnetic field 
orientation. The [110] Eixis of the sample was oriented perpendicular to the 
magnetic field, and changes in the magnetic field orientation were accomplished 
by rotating the sample about this [110] axis. In this way, the major crystalline 
symmetry axes ([001], [110], and [111]) could be aligned with the magnetic field. 

After a sample was studied, it was subjected to thermal damage, and each isotope 
was studied again in the same manner. 



In order to introduce a known defect concentration in the samples , 
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each was thermally damaged following the method of Potts and Pearson . In that 

work the crystals were held at a constant elevated temperature (1000*^C to 1200^C) 

in an evacuated vessel for 2^+ hours and then quenched to O^C in less than 0*2 

seconds. Vacancy concentrations were then determined by the Kossel line 

technique*^ ’ . In our experiment much lower prequench temperatures were used 

(500*^C to TOO^C), thus introducing lower defect concentrations. In addition, due 

to the size of our samples, a rapid quench would cause macroscopic damage, so a 

nnich longer quenched period of 10 to 15 minutes was used. Ionic mobilities in 

GaAs near room temperature are quite low , so the vacancy concentration caused 

by sublimation of As atoms is not expected to depend upon the period of the quench. 

However, Frenkel pairs may recombine quickly, and a longer quench, period might 

substantially reduce their concentration. Room temperature annealing of damaged 

samples was observed by Potts and Pearson, and it was interpreted, in part, as 

due to the recombination of Frenkel pairs. No annealing was observed in our 

samples, so any Frenkel pairs must have already recombined. In addition, after 

the quench a metallic film remained on the vycor vessel. Chemical analysis 

failed, but from vapor pressxire data we are convinced that it was As. It is 

concluded that this slow thermal quench generates As monovacancies, 

NMR measurements were made with a pulsed NMR system shown in block 

diagram form in Fig. 1. All measurements were made at a fixed frequency of l4.4 

MHz, with the appropriate magnetic field supplied by a Varian V-3900 electromagnet 

and power supply equipped with Mark I Fieldial stabilization. The rf system, 
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with few exceptions, is very similar to that described by Ellett,et al. as 




figure 1 Block diagram of pulsed NMR apparatus. 





spectrometer B. Oiir gated power amplifier, however, is an NMR Specialties Model 
P-103R. Also, the signal duplexer and preamplifier are those described by Lowe 
and Tarf,^^ modified for a 50 ohm system. A single coil, series tuned probe was 
used. At resonance the probe impedance was 50 ohms with a g.uality factor of 
approximately 5 • 

Since the sample is piezoelectric it produced an unwanted transient 
signal following the rf pulse. Without special precautions this distorted 
approximately the first 100 microseconds of the FID. It was suggested that such 
noise may be damped by placing the sample in contact with a fluid of appropriate 
viscosity. However, it was easier, and just as effective, to surround the 
sample by a Faraday shield. The mylar strip was removed from a metallized mylar 
capacitor, enclosed in insulating tape, and wound closely around the sample 
coaxially with the rf coil. This reduced both the amplitude and the duration of 
the piezoelectric noise by an order of magnitude. 

Control of the experiment timing is provided by a pulse programmer 
designed and constructed in this laboratory. Its operation is based upon a 
256 bit random access memory divided into 16 eight bit words. On execution, 
each word performs a single step of a pulse sequence, either a pulse or a delay. 
Pulse and delay durations can be controlled either digitally (using a 10 MHz crystal 
oscillator) or continuously (using 96OI monostable multivibrators). Two internal 
step counters as well as a minicomputer interface allow the programmer to control 
almost any pulsed NMR experiment. 

The signal to noise ratio (S/N) of the FID following each 90° pulse was 



approximately 10, too small for accurate determination of the linewidths. So 
5000 FID*s were signal averaged to yield a S/N of over 700. 

The FID undergoes analog-to-digital conversion in an 8 hit converter, 
a modified Computer Labs Model HS-802, at a rate of up to 5 MHz. Because this rate 
is too great for direct computer storage, the data for a single FID (1024 eight 
bit words) are stored in a buffer memory. Subsequent transfer to a Digital 
Equipment Corp. PDF 8/e minicomputer takes place during the NME longitudinal 
relaxation time, T^, at a rate set by the con^uter. The computer maintains 
storage of the most recent FID and the accumulated sum. Either 
of these may be displayed on an oscilloscope during T^. When 5000 FID's have 
occurred, the computer transfers the accumulated data onto 9 track, industry 
compatible magnetic tapes. The analysis of the data was performed with an IBM 
360/50 computer. 



IV . MALYSIS 


Lowe and Norberg^^ have shown that the moments of the NMR line can 
be determined from the time derivatives of the FID at t=0 by using the relation 








c-0" 


gco)V oievt'-o 


( 17 ) 


where is the 2n^^ moment and g(t) is the FID. However, the t 0 

point Is precisely in the center of the rf pulse^^^ PIP derivatives 

cannot be .ensured there directly. One method of obtaining second moments from 


the FID is to fit the visible part to an analytic function and to extrapolate 
back to t=0 in order to evaluate the derivatives. The accuracy of this method, 
of course, depends upon the ability of the function to describe the FID. 

We used a least squares fitting routine^^ to fit the FID with two 
functions. The Abragam function 

(18) 
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was found to approximate the F^^ free induction decay in CaF^ very well. Our 
data appear qualitatively similar to the fluorine FID. However, when the 
Abragam function was fit to GaAs FID's, there were systematic deviations. Fig. 2 
shows the generally unsatisfactory character of the fits of this function. 
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FID and liest functional fit using Abragast function. 


has spin ^ and so zero quadrupole moment- Q?hus one vould expect the 
GaAs FID shapes to be modified by (luadrupolar broadening* For half odd integral 
spins the first order effect of an efg is to leave the transition between the 
1/2 to “1/2 energy levels unchanged and to split the other transitions into 
satellites* If the efg is random this leaves the central component of the line . 
imbroadened while the remainder of the line broadens. Thus the FID has a fast 
decaying component arising from the satellites and a slower decaying component 
arising from the central component# In the case of spin 1=3/2 nuclei, the satellite 
lines make up 60 ^ and the central component lines the other k0% of the total 
resonant absorption. Thus the FID can better be described by a sum of two 
Abragam functions# The function 


describes the GaAs FID very well, as shown in Figs. 3 and U, The angular 
frequency second moment of a FID described by Eq# (19) is * 

Based on the above reasoning, more adequately discussed in the Appendix, 
one expects to find and^^S* However this is not always 

the case. We find that the ratio expected value with 

increasing quadrupole interaction# The quadrupole moments as shown in Table I 
otey the inequality In all samples the As'^ FID (for0 7^'n/2) is 

69 

described by Eq. (19) with Ij^/Ig) 2/3. This is also true of the Ga and, to a 


.71 


lesser extent ,the Ga FID in samples with large defect concentrations. 


One reason for this departure is our choice of the Abragam function as 


the functional form of the FID in the absence of the quadrupole interaction. 



TABLE I. RELEVANT NUCLEAR PROPERTIES®- 


Isotope 

Nat-oral 

Abundance 

{%) 

I 

Magnetic 

Moment 

(/»«) 

Electric 
Quadrupole 
Moment 
(barns ) 

^lU 

(10^° cm~^) 

Ga"^^ 

39.8 

3/2 

2.551^9 

.12" 

1.06® 

Ga^9 

60.2 

3/2 

2.0108 

.19'^ 

1.28® 

As'^5 

100. 

3/2 

1.431^9 

.29^ 

1.72® 


a. F. Bovey atid A. Tiers, NMR Tables (Wiley-Interscience , Inc. New York 196?), 

5th ed. 

b. The magnitude of coupling constant relating efg to applied electric field is 
adjusted for the given. values of quadrupole moment. 

c. G. F. Foster, Phys. Rev. 1^8 (1952). 

d. V. S. Korolkov, A. G. Makhanek, Opt. Spectry. USSR (English Transl.) 12,87 (1962). 

e. D. Gill and N. Bloembergen, Phys. Rev. 129, 2398 (I963). 











The Abragam function is successful in describing the FID of a spin 1/2 system 
broadened by the dipolar interaction only. The presence of the like nuclear 
exchange interaction can alter the shape of the FID substantially without an 
effect on the second moment . Further, in our lineshape calculation in the 
Appendix we do not allow the width or shape of the Abragam function to be 
affected by the guadrupolar interaction. Clearly the presence of a large efg 
can cause parts of both the like spin dipolar and like spin exchange interactions 
to become non-secular and may substantially alter the lineshape as well as the 
second moment. 

Certainly part of the deviation of the experimental lineshape from 
our simplified theory is due to o\ir inability to observe enough of the initial 
part of the FID so that the fitting routine may determine the parameters 

75 

accurately. From our theory, the satellite component of the As 

“3 

FID at 6 = 36 ^ in a sample of defect concentration 10 cm should have a 

half width of 23 ^sec • Thus we can detect only three quarters of the initial 

75 75 . ♦ . 

Intensity of the As FID, Clearly the As data are unreliable in this region 

of defect concentration for this magnetic field orientation, the [111] direction* 

In order to measure the broadening in samples of large defect concentration, second 

moment measurements were made near theG=1T/2 orientation where the initial part 

of the FID is slow enough to observe accurately. However, even here there are 

fitting difficulties in samples with large defect concentrations. 

As a result of the deviation of the fitting parameters of Eq. (l9) from 

their expected behavior, one can not always interpret these parameters physically 



as in the Appendix. However, "based on the agreement "between the FID and Eq. (19) 
the second moments are believed accurate. We use only the second moments in the 
following analysis. 

It was shown in Section II that a rotation of the magnetic field about the 
[iTO] direction results in the second moment varying as 

<aH^) = C+DfCs'). ( 20 ) 

The contributions to the isotropic part C are from the dipolar and exchange 
mechanisms, while the contributions to D are from the dipolar, pseudodipolar and 
quadrupolar mechanisms . Parts C and D were determined by fitting the 
measured second moments at known © to Eq.. (20). Generally, sixteen orientations 
covering a rotation of 120*^ were used to determine C and D. However, in samples 
with large defect concentrations only the orientations close to © =Tt/2 were 
used, as few as six orientations. Figs. 5> 6 and 7 show the second moments versus 
orientation for several samples . Table II contains the coefficients C and D for 
each isotope in each sample. 

Given these isotropic and angular dependent parts, the second moment 

contributions from the different broadening interactions may be identified. The 

IT 

pure van "Vleck dipolar contributions are well known , and they are for each 
isotope : 




(21a) 

= .as + uiice) 


(21b) 

- .a^+ 


(21c) 


where the subscript denotes the broadening mechanism, and the superscript the 

7 

isotope. NAR data on GaAs yield exchange coupling constants. Including only 
nearest and next nearest neighbor interactions, the unlike exchange contributions 



TABLE II, MEASURED SECOND MOMENTS AS A FUNCTION OF 
MAGNETIC FIELD ORIENTATION. THE SECOND MOMENT 

2 

IS EXPRESSED AS C + D f(0) IN UNITS OF GAUSS 


Sample^ 

Ga 

C 

D 

C D 

As'^^ 

C D 

MTU 

.58±.02 

•68±,o6 

. 03 

l.l8i.o8 

1.44±.07 

4.731.36 

M2U 

,6U±,03 

.631.07 

.641.03 

1.221.10 

1.321.07 

5.581.66 

M3U 

.65 ±.o 4 

.691.10 

.721.03 

1.301.11 

1.371.01 

4.631.04 

MTD500 

.581.03 

♦77 ±.o 8 

.631.01 

1.791.02 

1.401.10 

7.301.85 

M2D550 

.6^1.03 

1.301.10 

.66i.o4 

3.921.19 

1.231.08 

17.211.23 

M3D600 

.69±. 01 

1.841.03 

.721.04 

5.051.51 

1,151.08 

19.411. 60 

mtd 6 oo 

.631.04 

1.981,13 

.721.04 

6.081,45 

1.301.14 

24.911.61 

M5D700 

.651.04 

2.401.23 

.881.07 

7.461.58 

1.66i.l4 

21.211. l4 

EMC (Si) 

.721.05 

3.181.29 

.951.08 

7.341.73 

2,831.23 

25.112.34 


a) Sample identification M2D500 means a Monsanto crystal quenched from 500 C. 

Undamaged samples are labeled with a U. EMC(Si) is 'a Si doped sample^purchased 
from Electronic Materials Corp. with a carrier concentration of ^xlO cm ♦ 
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Figure 7 As ^ second moments as a function of magnetic field orientation for 


to the second moments are 



.e5 


(22a) 




(22t) 


.58 


(22c) 


Su'btracting "the contributions (2l) and (22) from our data leaves only 

the contributions from the quadrupolar and pseudodipolar interactions • Our 

analysis predicts that these contributions vill have no isotropic part, and, as 

can be seen from Table II and Eqs. (21) and (22), the isotropic remainder is 

negligibly small for the two Ga isotopes. In samples of large quadrupolar 

broadening, however, the isotropic remainder is larger than predicted on the basis 

of dipolar and exchange contributions. 

In order to separate the contributions from the quadrupolar and 

pseudodipolar interactions we use the fact that the pseudodipolar broadening in 
2 

gauss is the same for both Ga isotopes. This is because only the four nearest 
neighbors, all As nuclei, contribute significantly to the pseudodipolar broadening. 

Then the difference in the remaining angular parts of the Ga second moments must be due 
to quadrupolar broadening. Both Ga isotopes experience the same random efg, so 




respectively ; then 




(24) 


The relevant nuclear properties are shown in Table I. The values of 
^ik been adjusted for more recent vaJ_ues of electric quadrupole moment. 

Gill and Bloembergen give different numbers for the values of 

for the two Ga isotopes. Physically we see no grounds for this difference, and 
thus a weighted average value of 1.2xl0^^cm ^ is used for both. The 

69 69 

Ga value is weighted heavier because the experiment on Ga was probably more 

71 

accurate, due to the better S/N and greater quadrupole shift than Ga . Another 

ik 

experimental measurement of has been perfoimied . In this experiment 
was measured dynamically by observing the nuclear transitions caused by an rf 
electric field. The evaluation of R^j^ involves more manipulation than the 
straightforward experiment of Gill and Bloembergen. For this reason the static 


measurements of R^^^ are used here. 


Once and have been established we may find 

the Ga pseudo dipolar contribution 

<aH^)p = - (K -0 (kOJ - ^ (25) 


Comparing Eqs . (25) and ( 9)9 we obtain an a quadratic expression for the 

^ . 71 69 

pseudo dipolar coupling constant B in terms of D and D_ . There are two 

EiM R R 

solutions for CLnd. this experiment cannot distinguish between them- In the 

case of GaAs both solutions are negative in sign. 

The pseudodipolar interaction is proportional to the magnetic dipole 
moments of the interacting nuclei^. Thus we may write 




(26) 



where .U and .6 are the relative abundances of Ga"^^ and Ga^^, respectively. 


The 

T5 


quadrupolar second moment contribution is Just the remaining angular part, , 


less the pseudodipolar contribution . With the quadrupolar 

contributions to the second moments determined, we may calculate the concentration 
of charged defects. For nuclear spin I =* 3/2, Eq. (l6) reduces to 




(27) 


The defect density A is foimd from Eq. (2?) by identifying with 

the measured quadrupolar second moment. Table III contains the quadrupolar second 
moments 5 the pseudodipolar coupling constant, and the charged defect densities 
for each sample. 

Agreement of the defect densities as measured by G-a and As resonances 
is obtained only if we assume that the defects are preferentially located on the 
As sublattice. When this is true, the As nuclei are generally further from a 
defect than the Ga nuclei and will experience a smaller quadrupolar interaction. 

Thus the value of the lattice sum in Eqs. (l6) or (27) depends upon which sub- 
lattice contains the defects. The defect densities listed in Table III axe 
calculated assuming As monovacancies . 

Another assumption in the calculation of the defect densities of Table III 
is that the effective charge of a defect is the electronic charge, or e*=e. We 

assume a model in which the GaAs lattice is mostly covalently bonded, but retains 

' . 27 

some ionic character. The calculation of Harrison yields a net charge of 

+0.87e on the Ga site and a net charge of -0.87e on the As site, where the electron 



TABLE III. QUADRUPOLAR SECXDND MOMENTS, DEFECT DENSITIES, AND PSEUDODIPOLAR 
COUPLING CONSTANTS. DEFECT DENSITIES ASSUME AN ARSENIC MONOVACANCY 

OF CHARGE e. 


Sample 




Ga 


71 


Ga 


(gauss^) 

69 


As 


75 


Defect^- __ 
Density (10 ^cm 
As Measured By 
Ga As 


Pseudodipolar 

Coupling 

Constant 

%N 


M7U 

.171.03 

M2U 

.181.04 

M3U 

.201.05 

M7D500 

.341.03 

M2D550 

.861.07 

M3D600 

1.121.17 

M7D600 

1-331.15 

M5D700 

1,671.26 


.67±.13 3.85±.Hi 

♦ 73±.i 6 i+.75±.70 

• 8l±.20 3.83±.3i+ 

1.36±.11 6.62±.89 
3.it8±.28 i 6.6±1.28 

U.51±.68 i8.5±1.66 

5.361.62 23. 711.70 
6.7311.04 19.8±1.4U 


.21±.o4 .241.03 

.221.05 .301.05 

.201.06 .231.02 
.421.04 .421.06 

1.071.09 1.04±.08 

1.391.21 l.l6i.ll 

1.651.20 1.491.11 
2.071.32 1.251.10 


-1.54, -.401.05 
-1.52, -.48+. 06 
-1.53, -.471.08 
-1.47, -.531.06 
-1.48, -.521.09 

-1.57, -.431.08 
-1.65, -.351.12 
-1,70, -.301.18 


has charge -e ♦ If ve remove a neutral As atom from the lattice and allow 
electronic charge shifts on only the four nearest neighbors, then each neighboring 
Ga ion acq^uires a net charge of 0.55e, Now there will be a certain ionization 
energy for these four neighboring Ga ions to collectively give up one electron 
to the conduction band and in the process acq^uire an additional #25^ charge each. 
Then each Ga ion has a charge 0.9e> very nearly the 0.87e each had before the 
vacancy was produced. This small change in effective charge implies a small change 
in the bonding orbitals . On this basis we expect that the activation energy for 
ionization is small. The induced efg is sensitive to the change 
in the charge distribution. This difference is just 0-8Te localized at the 
site of the As vacancy and +0.03e at each neighboring Ga site. This appears 
very much like a point charge defect to all nuclei, even the nearest Ga neighbors. 



V. DISCUSSION 


a) Pseudodi polar Interaction 

The presence of the pseudodipolar interaction in III-V compounds has 

0 

been recognized only recently . One reason for this is that the nuclei in 

31 

these compounds generally experience quadrupole interactions ( P is the only 
exception). Since the magnetic field orientation dependence of the second 
moment contributions from each interaction has the form C+Df(0), the quadrupolar 
broadening can easily mask the effect of the pseudodipolar interaction. There 
would have been no clear evidence of the pseudodipolar interaction in this experiment 
had not our crystals been good enough so that the second moments of the Ga resonances 
were smaller than those calculated from the dipolar interaction alone. We now 
expect that the pseudodipolar interaction will be found in all III-V compounds. 

There have been several theoretical treatments of electron coupled 

4—6 

interactions, both from a band structure approach and using localized bonding 
13 26 29 

models. ^ In order to calculate the coupling constants from band theory, 

greatly simplified band structures were used. These calculations are very coarse. 

The localized bonding model calculations predict the negative sign of . observed 

1 O 

in TlCl InP , and now in GaAs. However, these calculations do not agree 


quantitatively with our data or the InP data. 


13 


In the treatment by Clough and Goldburg , it was noted that while 



the absolute magnitude of depend upon the degree of covalency A 

of the bond, the ratio B. ./K. .1 is independent of A . In their model they 

’ I ij ijl 

find that this, ratio has a maximum value of 1/2 for purely p orbitals. Larger 

g 

values have been measured in III-V compounds. Englesberg and Norberg observed 
a ratio of 1.8 in a single crystal of InP and 1.5 in a powder sample. 

We measure a value of .73 or 2.2 in a single crystals of GaAs, depending upon 
which solution of Eq. (9) is. physical. 

Clough and Goldburg developed their expressions for TlCl crystals. 

One should not expect that the same model can accurately describe III- 
V compounds. Their method begins by assuming a fixed degree of covalent 
character to the bonds. This covalent part of the bond is then represented by 
a variable admixture of s and p orbitals on each of the atoms contributing to the 
bond. They then find* using second order pertxxrbation theory, an expression for 
the Is. ./X. .1 ratio as a function of the amount of s and p character in the 

I ij ijl 

bond. 


The application of their model to GaAs has two major defects. It does 

not account for the d state admixtures in the wave functions that are known to 

be important. The second problem is that in polar semiconductors the wave 

27 

function of valence electrons moves toward the anion (As)» while the wave function 
for the conduction band shifts toward the cation^^ (Ga)* Since the exchange and 
pseudodipolar interactions are related to the overlap of the valence and conduction 
band wave functions, a model like Clough and Goldburg’ s, that assumes a fixed 
covalent character to the bond for both the valence and conduction bands, will miss 
the consequences • of the wave function shift between, the cation and the anion. A 
modification of their calculation that accounts for the relative displacements of 
valence and conduction electron wave fxmctions and the d state admixtures should 



be more accxirate . 


b.) Quadrupolar Effects 

An independent method of estimating the defect concentration is 

desirable, and the experimental data of Potts and Pearson enable one such 

determination. Their minimum quenching temperature was 1000°C, generating a 

defect density of 2 xlO"^ cm . Thus we must extrapolate over three orders 

of magnitude to compare with the defect density expected following a quench from 

a temperature in the 500°C to T00°C range. In Fig. 8 the log of the defect 

density is plotted against the reciprocal quenching temperature, T The measure- 

A. 

ments of Potts and Pearson are shown, as well as the dashed line showing the 
extrapolation based on their entropy and enthalpy values for the formation of 
a vacancy. 

It is clear that our observed concentrations are greater than those pro- 
jected. However, Potts and Pearson expressed some uncertainty in the proportionality 
constant relating the rms deviation in lattice parameter and the concentration of 
monovacancies • This puts the accuracy of their measured vacancy concentrations in 
doubt by a multiplicative constant. Further, the reliability of an extrapolation over 

three orders of magnitude is questionable. In Fig. 8 it can be seen that a small 

19 -3 

error in the activation energy found at a defect concentration of 10 cm makes a 

15 -3 

substantial difference when extrapolated to 10 cm . In addition, there is no 
a priori reason to believe that the same physical mechanism is dominant over the 
entire temperature range. 

We believe that the As vacancies are caused by the As atoms lost to 
the equilibrium vapor phase dixring the heating process. The equilibrium As^ and 
As^ vapor pressures are known^^ as a function of temperature in the 6T5^C to 925*^0 
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Also shown are the defect measurements of Potts and Pearson. The 
dashed line Is an extrapolation of their measurements to smaller 
defect densities. The solid line shows our calculation of defect 


densities using vapor pressure data* The error bars shown are typical 
of that temperature. 



range. Thus we can calculate the number of As atoms removed from a sample heated 

in an initially evacuated tube of a given size. The ratio of this number to the 

sample volume is the concentration of As monovacancies produced by this process. 

This calculation is shown in Fig. 8 as the solid line. Using the condition 

that each defect of charge e is located on the As sublattice, and using the low 

frequency dielectric constant for €. , our observed defect concentrations are 

smaller than the vapor pressure calculation. This calculation is meant for order 

of magnitude comparison, and it assumes that the ratio of tube to sample volumes 

was 10 for all samples. Generally, our tube volumes were larger, so the solid 
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line in Fig. 8 should represent the minimum number of defects produced. 

If the sample behaves as an infinite source of As vapor, supplying 
whatever is necessary to reach equilibrixun at a given temperature, then the defect 
density present following a quench depends additively upon that prior to the 
quench. In order to determine the temperature dependence the defect density in 
each sample prior to the quench must be subtracted from the measured density 
following the quench. Our data in Fig. 8 have been adjusted to account for 
prequench defect densities. 

The slope of our lowest temperature data is in agreement, within 
experimental error, with the vapor pressure activation energy. Using the two 
most reliable temperatures, 500^0 and 550 ^ 0 , we measure an activation energy of 
1.5 ± 0.1 eV. The activation energy reported by Potts and Pearson is 2.0 eV, and 
the vapor pressure data yield 1.4 eV. 

There are severaJL possible explanations for the fact that our measure- 
ments of defect concentrations are smaller than those calculated from vapor 



pressure data* One is the dependence on the experimental value of the product 
2 2 

^lU^ * small error in the values of or Q can maite a substantial 

error in our analysis of the defect densities. The error bars shown in Fig. 

8 do not account for errors in these quantities. 

Another explanation is that the thermal damage did not produce as many 

vacancies as predicted. This could have occurred in two ways: l) the sample 

was not in thermal equilibrium with the vapor, so the vapor phase did not consist 

of as many As molecules as predicted; or 2) the vapor pressure data are wrong. 

19 

Using self diffusion data it is estimated that the sample will reach thermal 
equilibrium in approximately 10 hours, and also it is unlikely that the As vapor 

can re-enter the lattice during the quench. Further, the vapor pressure data 

31 

accuracy is well established by several workers, so we have confidence in the 

number of As atoms in the equilibrium vapor phase, it is concluded that 

at least as many vacancies were produced as predicted hy the vapor pressure data. 

This experiment measures the charged defect concentration. The effect 
of the strain caused by an uncharged impurity is not sufficiently strong to be 
noticed when charged defects are also present. Another possible 
explanation of our apparently low defect concentrations is that not every 
As vacancy acquires a charge of magnitude e. This can happen in two 
ways. One is that each As vacancy assumes an effective charge e* such that e*<e. 
However, we have already shown that the model in which e*=^e discussed in Section IV 
is reasonable. The other is that the activation energy for- ionization of the defect 
is such that at room .temperature not all the impurities are charged. An experiment 
investigating the temperature dependence of the NMR second moment is expected 
to check this latter possibility. 



Our analysis uses the induced quadrupolar interaction, and the efg 

_2 

associated with a charged defect decreases as r because we have assumed that 

the defect appears as a point charge* However, electronic screening or defect 

-2 

complexes may exist, causing the efg to decrease faster than r . We have 
calculated the Debye-Huckel screening length for the Monsanto sample and find it 
to be negligible, so screening by free electrons is not important- Defect 
complexes can, however, make a substantial effect- For example, an As vacancy 
may pair with an acceptor to produce a localized electric dipole moment - 
In this case the efg will decrease as r the effective broadening range of the 
defect will be reduced, and many more impurities will be necessary to produce a 
given c^uadrupole contribution to the second moment* A double resonance experiment 
detecting nearest nuclei of the impurity as the rare spins will address the 
question of defect complexes. 

c) Experimental Limitation 

It is possible that the low values of defect concentration measured in 

samples quenched at 600^C and TOO^C are accurate. This effect could be caused by 

vacancy clustering or vacancy -acceptor pairing. However, it is more likely that 

that these low values are due to the current limitations in our experiment and 

analysis. There are several trends in the data that indicate that the present 

15 -3 

measurement of defect density errs if concentrations are larger than ^10 cm 

One trend, discussed in Section IV, is the observed deviation of the FID functional 

fitting parameters in Eq. (19) from physically interpretable quantities. This 

75 69 Tl 

deviation is always present in the As and Ga data. It appears in the Ga data 



for samples quenched from 600^C and higher- We also find that the increases 
in the second moment caused hy thermal damage at 600 ^C and TOO^C are proportionately 
smaller for nuclei with larger quadrupole moments- That is, the ratios of the 
increase in the second moment to the increase in the Ga^^ and second 

moments are greater than expected. These ratios should he fixed by the nuclear 
properties. It is clear that the fast component of the FID decays so quickly in 
the presence of these large quadrupole interactions that our fitting analysis is 
not sensitive enough to it. With our analysis this results in a deviation in the 
pseudodipolar coupling constant and a disagreement between the defect densities 
as measured by Ga and As resonances - 

Additional evidence that this experiment cannot measure relatively large 

69 

defect concentrations is seen in the data for the Si doped sample- The Ga 

second moments are much smaller than one expects from the large second moments of 

?1 75 

Ga . To an even greater extent, the same is true of As . Assuming that the 

pseudodipolar interaction is the same in all GaAs crystals, and assuming that the 

71 

ionized Si donors are located on Ga sites, then the Ga data indicate a charged 

15 — 

defect density of 6.1+xlO cm , where we have accounted for the Debye-Hiickel 

screening. Hall effect measurements yield a carrier concentration of 4x10 cm 

so there are at least that many ionized impurities. With our current experimental 

limitations, this technique is too sensitive to accurately measure charged defect 

15 

concentrations greater than 10 cm . 

We conclude that our present experiment does not dependably measure 

2 

frequency second moments larger than 5 kHz - This corresponds to a FID half width 
of approximately 70 ;isec. In the Ga^^, Ga^^, and resonances this corresponds 



to charged defect densities of 2.5x10^^ cra”^, 9x10^^ cm~^, and 2x10^^ cm"^, 
respectively. In order to make equivalent cw measurements, excellent S/N is 
required at least as far as 5 kHz from the center frequency. With improved 
damping of the piezoelectric noise and with spin locking techniques we expect 
these upper limits to he extended to charged defect densities near 10^"^ cm"^. At 
these concentrations we expect that second order quadrupolar effects will he 
important, and that the lineshape will have changed substantially. 

The lower limit of measurable charged defect concentration would occur 

where the quadrupole contribution to the As second moment is just measureable , 

say 0.1 gauss . This corresponds to a defect density of 6x10 cm on the As 
12 “3 

sublattice or 3x10 cm on the Ga sublattice, representing a defect to atom 
ratio of the order of 10 . When this method is applicable , it rivals the 

sensitivity of neutron activation measurements, and in addition, the WMR method 
can observe vacancies. 

In summary, while questions remain about the effective charge and the 
volume dependence of the defect densities, it is clear that the NMR lines in GaAs 
are extremely sensitive to impurity concentrations. In fact, they are so sensitive 

if) -3 

that in samples of carrier concentration of 10 cm“^ this first order line shape 
and second moment analysis is insufficient to treat the problem. Using this 
undoped semi -insulating single crystal, we have observed a negative coupling 
constant for the pseudodipolar interaction and developed a method by which the 
pseudodipolar, dipolar, and quadrupolar second moment contributions can be 
identified. We have made two estimates of the defect concentrations produced by 
thermal quenching, using independent sets of e3q)erimental data. Our measurements 



lie between these two calculations, above the Potts and Pearson projection and 
below the vapor pressure estimate. That our method of analysis is correct is 
adequately demonstrated by two facts. One is that the pseudodipolar interaction, 
within experimental error, remains independent of the total second moment, i.e, 
thermal damage. The second is that the charged defect densities as measured by 
Ga and by As resonances agree for samples with small concentrations. 
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APPENDIX 


The objective of this appendix is to compute the WMR line shape due to 
the quadrupolar interaction using a continuiun crystal model. Prom Eqs. (3g), (12) 
and (13) the satellite frequency shift of a spin 3/2 nucleus located a distance ^ 
from a charged impurity is 


aV = 


-oc • r 



(Al) 


\diere 


ot ~ 


€ K 




(A2) 


We calcrilate the volume XlL about a single impurity containing nuclei with frequency 
shifts greater than some aV* 

The volume il is bounded by a revolution about the ot direction of a curve 

<vw 

described by 

r® Z - c?^Cos<}>/4V ^ (A3) 

where (f> is the angle between r and ^ . For ^ =0 there is a minimrun frequency shift 
for the nuclei located on the crystal surface, a distance r^ from 
the impurity* The maximuin frectuency shift is 5 f'or the nearest 

neighbor nuclei located at a distance ^ from the impurity* The volume within the 



nearest neighbor sites and outside the sample are excluded from . These 
exclusions eliminate unphysical infinities. The line shape will he shown to he 
independent of hut not of a. Thus the sample dimension has no effect upon 
the line shape. Further, the near neighbor distance a may he used to determine the 
location of the defect sites. The volume XL is found by a simple integration, 
and the result is 


XlM 


( I 5 — O, ^1“ ■§ Av^ 

L 0 , Av>„<AV> 


The number of nuclei M 0i’\>) shifted by an amount greater than A'v* is 

5 

(3AyH-n<^)£}(AV), where N is the number of resonant nuclei in the sample. Then 
the distribution g(AV) of frequency shifts caused by a single impurity is 



dNs 


(A5) 


Combining Eqs . (a 4) and (A5) the expression for g(^V) becomes 


f \ - A\><AV, 




- 3a/ 


(AVA^>)^-(^VA'^'^A^l^A^>:<A^>.„ (A6) 


//a 


1^0, < 


AN> 



Next the frequency shift distribution due to a number of impurities 

randomly distributed throughout the sample must be calculated. Due to the long 

range natvire of the induced efg, the satellite transitions of each nucleus are 

significantly affected by many defects, even at very low defect concentrations 
12 -^3 

( 10 cm ), Thus one may use the central limits theorem to calculate the frequency 
shift distribution due to a number of randomly distributed defects. According 
to that theorem, the distribution vill be gaussian with a second moment given by 


< = Ni ^ 

'Where is the second moment of g(AV). Using the relation a/r « 1, 

one finds 




(AT) 


OO 

where ^ is defined in Eq. (7), and is the charged defect density. 

Finally, the relation between the distribution of quadrupole frequency 
shifts G (basically an inhomogeneous broadening) , and the WMR line shape 

G(V) must be established. The quadrupole frequency shift distribution is a gaussian 
function with second moment 


For the spin 3/2 case, the magnetic resonance line shape is found 


through the expression 


G Cv)^ A\>)] (A9) 

_'io ‘ ' ^ 

\Aiere G^^^ 3 / 2 ^^^ shape function of the 3/2-1/2 transition caused 

by other broadening mechanisms (as modified by the quadrupole shifts). The 



is the center frequency 


first argument in the functional dependence of ^^^>2 3/2 
for the 1/2-3/2 transition* and the second is intended to represent the effect the 
shift may cause on the line shape. The shape due to other interactions can be 
modified by the shifts, e.g. terms in the dipolar hamiltonian that are secular in 
the absence of shifts will become non-secular when the shifts are large enough, thus 
modifying the width of the dipolar broadened line. 1/2 *^-3/2 -1/2 

corresponding shape functions for those transitions. The factors 0.3, 0.4, and 0.3 
are weighting factors that arise from the transition matrix elements for the three 
transitions. 

For the case of low defect concentrations treated here, we ignore the 
effect of the frequency shift upon the line shape. Further we assume that 

lineshape in the absence of the quadrupole interaction. Then Sq. (A9), in view 
of these approximations, becomes 


Gu) - O.mGoCV) tO.<b Jd(AV)C^Cv>+4v') G^CAv^) 


(AlO) 


^00 


If Gq(v^) is a gaussian with second moment then 


-i 

Gcv") = o.6[aTT(<dV^^+div’‘>,)] ‘ v>%)] t 

exp[- ^ 


(All) 


where is centered at the Larmor frequency. Thus the second moment <AV^> of 
the composite line is 


<AV‘> = +0.0; 


<4. 


(A12) 



The Fourier transform of G(\>) is the FID function. Thus the quadru- 


polar broadened FID should have the fonn of Eq. ( 19 ) with ~ 
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For comparison* we can examine Simdfors*5 theoretical model. He sets 
£ = [ ^ 6. - ^ then computes Eq_. (A3) hy setting r equal to 

^ 1/3 

one -half the average distance between impurities* r=l/2( 3/^’7T^x ^ » 

and averaging over ^ . This yields: 

<AV">, - (‘iw/’./sf k’Hce). 

Using his experimental values of he calculated values of that 

differ from other experimental values obtained by more direct means^^ by a 
factor of three. While Sundfors^s expression for has the correct 

angular dependence, it will not produce the correct numerical results. 



